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10.5 IF ESTIMATION OF FM SIGNALS IN MULTIPLICATIVE 
NOISE 0 

10.5.1 Random Amplitude Modulation 

Most IF estimation techniques, such as those presented in the previous articles of 
this chapter, assume that the signal of interest has a constant amplitude. While 
this is a valid assumption in a wide range of scenarios, there are several important 
applications in which this assumption does not hold. Indeed, in many situations 
the signal may be subjected to a random amplitude modulation which behaves as 
multiplicative noise. Examples include fading in wireless communications [1], fluctu- 
ating targets in radar [2] , and structural vibration of a spacecraft during launch and 
atmospheric turbulence [3]. In this article, we focus on non-parametric methods. 
In particular, we show that the Wigner-Ville distribution (defined in Section 2.1.4) 
is able to display the IF of a signal affected by multiplicative noise, and that this 
representation is optimal in the sense of maximum energy concentration for a lin- 
ear FM signal. For higher-order polynomial FM signals, the use of the polynomial 
Wigner-Ville distribution (PWVD), presented in Article 5.4, is shown to give opti- 
mal representations. Statistical performance of each case will be presented here. 

10.5.2 Linear FM Signal 

In this section, we study the case of a linear FM signal and assume that the multi- 
plicative noise is a real-valued process. 

10.5.2.1 Optimality of the Wigner-Ville Spectrum 

First we show that the Wigner-Ville spectrum (WVS) is optimal, in the sense of 
IF localization, for the time-frequency analysis of linear FM signals affected by 
multiplicative noise. Consider the signal y(t) given by 

y(t) = a(t) ■ z(t ) (10.5.1) 

where a(t) is a non-zero-mean real-valued stationary noise and z(t) is a deterministic 
FM signal given by z(t) = exp {j<f>(t)j. For a linear FM signal, </>(t ) is a second- 
order polynomial. Using the expectation operator notation, the autocorrelation of 
the signal above can be expressed as 

ICy(t,T) = E \y(t - |) y*(t + §)] 

= E[a(i-§)a(t+§)] .{ Z {t-l)z*{t+l)} 

= U a (T)K z (t,T). (10.5.2) 
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The WVS of y(t), which is defined as the Fourier transform of K y {t,r ) [4], can be 
expressed as 

Wy(t,f) = ^f{ICy(t,T)} (10.5.3) 

= Tr^j {z{t - i) Z*(t + §)} *f Tr^f {E [a(t - i) a(t + §)] } (10.5.4) 

= W z (t,f)* f S a (f) (10.5.5) 

where *f is the convolution operation in the frequency space. 

If we express the non-zero- mean random process a(t) as a{t) = fi a + cio{t), where 
/ i a is a constant mean of a(t) and ao (t) is a zero- mean noise with autocorrelation 
i? ao (r), we can re-write the WVS of y(t) as 

Wy(t, f) = nlw x {t, f ) + S ao (f) *f W z (t , /). (10.5.6) 

For the case of a linear FM signal, the Wigner-Ville distribution (WVD) is given in 
Article 2.1 by [4] 

W z (t,f) = 6(f-m ), (10.5.7) 

where /)(f) is the signal IF and (5 is the Dirac delta function. In this case, we obtain 

Wy(t, f ) = n 2 a S{f - fi(t )) + S ao (f - Mt)). (10.5.8) 

Note that Eq. (10.5.8) exhibits the presence of a spectral line at the frequency /;(t) 
for all time instants. This means that, theoretically, the WVS always localizes the 
IF of a linear FM signal. This makes it a powerful tool in the analysis of linear FM 
signals affected by multiplicative noise. Also note that when y a = 0, the WVS will 
not exhibit a peak at the signal IF, indicating a breakdown of the WVS to analyze 
the noisy signal. 

As an illustration, consider a unit-modulus linear FM signal sampled at 1 Hz, 
whose frequency range lies between 0.1 Hz and 0.4 Hz. The signal length (in samples) 
is chosen as N = 511. This signal is multiplied by a real-valued i.i.d. Gaussian 
noise with a standard deviation a a = 1 and a mean equal to 0 and 1, respectively. 
Fig. 10.5.1 displays the WVS (one realization) of the noisy signal for both values of 
the mean. As expected, the WVS for the zero-mean case cannot reveal the signal 
IF; however, it can do so for the other case. In this last case, the peak of the WVS 
can be used to estimate the IF of the signal. In what follows, we will evaluate the 
statistical performance of such an estimator. 

10.5.2.2 Statistical Performance Evaluation 

Here, for a more complete study, we consider the presence of additive noise as well 
as the multiplicative noise. The objective is to derive the asymptotic variance of 
the IF estimator, based on the peak of the WVS, for this case. 

Let the discrete-time version of the noisy signal be 

y{n) = a{n) + w{n) , n = 0, . . . , N — 1. (10.5.9) 
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Fig. 10.5.1: The Wigner-Ville spectrum of a linear FM signal affected by a real-valued Gaussian 
multiplicative noise. The noise variance is 1 (both plots) and its mean is 0 (left plot) and 1 (right plot). 



The process a(ri) is considered to be a real-valued stationary Gaussian noise with 
mean and variance given by fi a and respectively. The complex zero- mean ad- 
ditive process w(n) is assumed to be stationary, white, circular and Gaussian with 
variance equal to crjj,. In addition, both noises are assumed to be independent. 

The WVS used to estimate the signal IF is defined, in the discrete-time domain, 
using the expression given in Article 6.1, as [4] 



Wz(n, /) = E 



M 

2 ^ y(n + m) ■ y*(n — m) e - - 7471 ^ 1 ” . 

_ m=—M 



(10.5.10) 



Using straightforward derivations, we can show that, for increasing window 
length (2 M + 1), the WVS converges in probability to [5]. We can 

also show that the IF estimator asymptotic variance is approximately equal to [5] 

= (2n)' 2 S w (2M + l ) 3 ( 2 + * + i) I 10 ' 5 ' 11 ) 

where S a = Ma/ a a anc I = Ha/&w- Note that: 

(i) When /i a = 0, the variance goes to infinity indicating that the WVS based 
estimator breaks down. This result confirms the analysis presented earlier. 



(ii) When ji„ = A where A is a constant and cr a = 0, i.e., the signal under 
consideration is just a constant amplitude linear FM signal embedded in noise, 
the variance expression can be-rewritten as 



Var (/i(n)) 



3 < 4[2 A 2 + cr 2 w ] 
(2tt) 2 A 4 (2M + 1) 3 ' 
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Fig. 10.5.2: Theoretical (dashed curve) and estimated ('+’) variances of a linear FM signal corrupted 
by multiplicative and additive noise. Left plot is for S w = OdB and right plot is for S w = 5dB. 



If, in addition, we assume high signal-to-noise ratio (SNR) (i.e., A 2 o 
the asymptotic variance expression reduces to 



Var(/i(n)) 



6(7 



2 

w 



(2tt) 2 (A 2 )(2M + l) 3 



which is similar to the result obtained in [6] . 

The above theoretical results were confirmed by Monte-Carlo simulations. 
Specifically, we estimate the IF estimator variance using 5000 realizations of the 
signal given by (10.5.9). In Fig. 10.5.2, we display the theoretical (dashed curve) 
and the estimated (“+’) variances plotted against S a for S w = 0 and 5dB, respec- 
tively. 

If the signal under consideration is not a linear FM but a higher-order polynomial 
FM signal, the WVS becomes inappropriate because it introduces some artifacts 
which might hide the real features of the signal and its peak based IF estimator 
is biased for such signals [6]. In this situation, a different tool is needed for the 
analysis. This is the topic of the next section. 



10.5.3 Polynomial FM Signals 

Compared to the previous section, the extension here is two-fold: (i) the signal 
considered is assumed to be a polynomial FM signal of arbitrary order and (ii) the 
multiplicative noise is no longer limited to a real-valued process but is assumed to 
be a non-zero complex circular Gaussian process. Based on this, the noisy signal 
y(t) is now written as 

y(t) = a(t) ■ z(t) + w(t) (10.5.12) 

where the stationary processes a(t ) and w(t) are both assumed circular, complex, 
Gaussian and independent with means and variances given by (/z 0 , a 2 ) and (0, a 2 , ) , 
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respectively. The noiseless polynomial FM signal z(t) is given by 

z(t) = = exp |j (10.5.13) 

where the a, are real coefficients and P is the order of the polynomial phase. Note 
that in the derivation below we do not require the knowledge of the coefficients a,; 
we only assume a polynomial FM signal. 

The IF of the signal z(t) is given in Article 1.3 by 



/*(*) 



1 d<j){t) 
27 r dt 



p 
i — 1 



and our primary objective here is to estimate fi(t) from the noisy signal y(t). For 
that purpose, we use the polynomial Wigner-Ville distribution (PWVD) defined in 
Article 5.4 as 







(10.5.14) 

(10.5.15) 



where q is an even integer which indicates the order of non-linearity of the PWVD. 
The coefficients c, and c_i (* = 1,2,..., q/2) are calculated so that the PWVD is 
real and equal to 

WP(n,f)=8(f-fi(t)), 

for signals given by Eq. (10.5.13). Note that the realness of the PWVD implies 
that Ci = — c_j. Also note that the WVD is a member of the PWVDs class with 
parameters q = 2 and ci = — c_i = 0.5. Full details of the design procedure may 
be found in [7] and Articles 5.4 and 5.5. 

The choice of the PWVD steins from the fact that it yields a continuum of delta 
functions around the IF for a given polynomial FM signal. This property implies 
that the peak of the PWVD can be used as an IF estimator for polynomial FM 
in a noisy environment. In [8], the statistical performance of this estimator was 
evaluated for noisy signals described by Eq. (10.5.12). It shows that this estimator 
is unbiased and its asymptotic variance is approximately equal to [8] 



Var(/i(n)) 



6(^ + <rg,)E?^fc? 

(27t)V 0 | 2 (2M + 1)3‘ 



(10.5.16) 



In the above expression, ri\ represents the number of the different coefficients c t in 
the PWVD kernel, while k t (for i = 1, . ... , m) represents the multiplicity of each of 
these coefficients cq. and (2 M + 1) is the window length considered in the PWVD 
discrete-time implementation. Note that: 
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(i) When y a = 0, the variance goes to infinity indicating that the PWVD based 
estimator breaks down. 



(ii) When /. i a = A where A is a constant and cr a = 0, i.e. , the signal under 
consideration is just a constant amplitude polynomial FM signal embedded in 
complex Gaussian noise, the variance expression can be-rewritten as 



Var(/0 



(2n) 2 A 2 (2M + l) 3 ' 



The above expression is exactly similar to the result obtained in [6], which 
treats constant amplitude polynomial FM signals only. 



10.5.3.1 Monte-Carlo Simulations 

To confirm the validity of the above theoretical results, we consider the IF estimation 
of a quadratic FM signal at the middle of the signal interval. The peak of the sixth- 
order PWVD, whose signal kernel is given by [7] 

K^(t,r) = [y(t + 0.62r) y*(t — 0.62r)] • [y(t + 0.75r) y*(t — 0.75r)] 
x [y(t — 0.87t) y*(t + 0.87r)], 

is used here as the IF estimator. The noisy signal y(t) is generated as suggested 
by Eq. (10.5.12). For this example, we choose the sampling period equal to T = 1, 
the signal length N equal to 129, the window length equal to (2 M + 1 = N = 129) 
and the noise variances to be equal (i.e., ay = <r, 2 ). In the simulations, the overall 
signal-to-noise ratio, defined as 

SNR W1 = 101og 10 (|^ a | 2 /(fj2 + <r 2 )), 

is varied in a 1 dB step from 0 to 15 dB. Monte-Carlo simulations for 1000 realiza- 
tions are run for each value of SNR„ ;i . The results of two different experiments, 
one performed for \y a \ = 0.01 and the other performed for \y a \ = 1, are displayed 
in Fig. 10.5.3 (left plot). We observe that, above a certain threshold, the estimated 
variances represented by ‘+’ (for \y a \ = 0.01) and ‘o’ (for \y a \ = 1) are in total 
agreement with the theoretical ones given by Eq. (10.5.16) and represented by the 
continuous lines (superimposed). 

Simulations run under the same noise conditions for other polynomial FM signals 
using the appropriate PWVD order also confirm the theoretical results presented 
above. One such case is when y{t) is a linear FM signal and the PWVD considered 
is the second-order PWVD. The results of this experiment are displayed in the right 
plot of Fig. 10.5.3. 



10.5.4 Time-Varying Higher-Order Spectra 



Time- Varying Higher-Order Spectra (TV-HOS) based on the polynomial Wigner- 
Ville distribution are defined as the expected value of the PWVD [9], namely 






(10.5.17) 
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Fig. 10.5.3: Results of the experiments for a quadratic (left plot) and a linear (right plot) FM 
signal corrupted by complex Gaussian multiplicative and additive noise processes. The continuous 
lines (superimposed) represent the theoretical variances while '+' and ‘o’ correspond to the estimated 
variances for \p a \ = 0.01 and \p a \ = 1. respectively. 



where Wz q \t,f) is the (/ th -order PWVD defined by (10.5.14). Interchanging the 
expectation operator with the integration in the PWVD, one obtains 







(10.5.18) 

(10.5.19) 



where K ( t , r) represents a slice of a time- varying g th -order moment function [9] . 
If the quantity ld 9 \t, r) is absolutely integrable, then, wi q \t, /) can be interpreted 
as a form of the time- varying higher-order moment spectrum. In [9], the authors 
showed that TV-HOS combine the advantages of classical time-frequency analysis 
with the benefits of higher-order spectra. To avoid the problem of non-superposition 
of the higher-order moments, the authors in [9] used higher-order cumulants instead. 

It is important to note that since non-stationary random signals are non-ergodic, 
the ensemble averaging above cannot be replaced by time averaging. In this situa- 
tion, local ergodicity has to be assumed. 

Readers interested in TV-HOS are referred to [9] for more details, including 
examples of the efficacy of TV-HOS in the analysis of random FM signals affected 
by multiplicative noise. (See also Section 14.5.4 and the references in [9].) 



10.5.5 Summary and Conclusions 

The Wigner-Ville spectrum (WVS) and polynomial Wigner-Ville distributions 
(PWVDs) are considered for the analysis of polynomial FM signals corrupted by 
multiplicative and additive noise. In the noisy linear FM case, the WVS is shown 






464 



Chapter 10: Instantaneous Frequency Estimation and Localization 



to give optimal IF localization. Accordingly, the peak of the WVS is proposed as an 
IF estimator. A statistical performance test shows that this estimator is very accu- 
rate even at low signal-to-noise ratio values. For the case of the noisy higher-order 
polynomial FM signal, the peak of the PWVD is shown to be a very consistent and 
accurate IF estimator. 
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